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EREFACE 

This  Msmorandiun  is  a  product  of  research  heing  conducted  at  RAND 
on  the  design  and  operation  of  simulation  models  for  studying  manage¬ 
ment  policies  and  other  problems  that  involve  con5)lex  systems  of  random 
variables.  Simulation  techniques  are  being  used  increasingly  in  manage¬ 
ment  studies,  for  instance,  in  examining  supply  and  maintenance  policies; 
and  it  is  desirable  to  make  the  technique  more  rigorous . 

The  text  is  a  companion  piece  to  M.  A,  Geisler,  The  Sizes  of 
Simulation  Samples  Required  to  Compute  Certain  Inventory  Characteristics 
with  Stated  Precision  and  Confidence,  The  RAM)  Corporation,  RM-3242-H?, 
October,  I962.  Special  statistical  methods  were  used  in  that  study  to 
con^te  the  san^le  sizes  for  specified  inventory  models.  In  this  study, 
the  methods  were  tested  by  applying  them  to  particular  inventory  cases, 
and  determining  how  well  the  actual  precision  and  confidence  obtained 
in  the  estimates  agreed  with  expectation.  Another  related  publication 
is  R.  W,  Conway,  Some  Tactical  Problems  in  Sigulation  Method,  The  RAND 
Corporation,  RM-52U4-HI,  October,  1962. 

Information  compiled  in  this  Memorandum  should  be  of  interest  to 
Air  Force  statisticians  and  research  people  engaged  in  applying  simu¬ 
lation  techniques  to  management  and  other  planning  problems. 

This  study  will  be  presented  at  the  Tenth  International  Meeting 
of  the  Institute  of  Meuiagement  Sciences,  Tokyo,  Japan,  in  August,  1965. 
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SUMMARY 

RAND  Memorandum  Ji2k2*  presented  l6  tables  of  sample  sizes  to  be 
used  for  estimating  certain  statistical  characteristics  (mean  shortages 
and  overages  per  time  period)  of  particular  inventory  models  by  com¬ 
puter  simulation.  The  tables  cover  a  vide  range  of  inventory  policy 
parameters;  the  inventory  models  used  include  the  zero  procurement 
lead-time  case,  plus  the  non-zero  procurement  cases  of  2-,  5“>  an<3. 
lO-period  lead-tiiass  with  exponential  demand.  To  compute  these  tables, 
we  had  to  make  several  important  assumptions  that  could  affect  the 
accuracy  realised  in  simulations  based  on  these  sample  sizes,  compared 
with  that  expected  from  statistical  theory.  We  undertook  this  current 
study  to  determine  how  well  realized  and  expected  accuracy  might  agree . 

The  test  consisted  of  applying,  to  1000  independent  samples,  the 
sampling  procedure  described  in  RM-5242,  using  several  inventory  policies 
and  the  four  procurement  lead-time  cases  mentioned.  (Even  thou^  we 
could  solve  the  zero  lead-time  case  completely  by  mathematical  analysis, 
we  examined  it  to  obtain  an  absolute  test  of  the  statistical  method. 

In  addition,  the  non-zero  procurement  lead-time  cases  were  tested,  using 
empirically  derived  standards . )  In  all,  we  tried  25  separate  inTOntory 
cases,  each  with  1000  samples  computed.  For  each  sample  size,  we  com- 
I*ited  the  mean  number  of  shortages  and  overages  per  period.  Then,  from 
the  1000  samples  for  each  inventory  situation,  we  could  compute  the 
percentage  of  samples  falling  within  the  specified  precision  and  confidence 
range  predicted  in  the  previous  Msraorandum. 

M.  A.  Geisler,  The  Sizes  of  Simulation  Samples  Required  to  Compute 
Certain  Inventory  Characteristics  with  Stated  Precision  and  Confidence, 

Bfce  RAND  Corporation,  RM-5214-2-ER,  October,  19^2. 


We  found  that  the  actual  precision  and  confidence  obtained  for 
each  of  the  inventory  policies  and  procurement  lead-time  cases  tested 
did  correspond  closely  with  expectation.  Thus,  it  was  expected  that 
95  per  cent  of  each  of  the  zero  procurement  cases  tried  would  be  within 
100  per  cent  of  the  true  mean  value.  In  the  shortage  calculations,  the 
values  ranged  from  9^*3  "to  97-6  per  cent,  and  for  the  overages,  the  - 
range  was  from  9?. 6  to  99*9  per  cent.  From  this  test,  we  concluded 
that  the  statistical  procedure  given  in  RM-32lt-2  is  valid,  and  produces 
reliable  estimates  of  mean  shortages  and  overages  per  period. 
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I.  INTRODUCTION 

The  availability  of  electronic  computers  has  made  it  possible  to 
extend  mathematical  and  statistical  results  beyond  those  attainable 
solely  by  analytical  techniques.  The  combined  use  of  analytical  and 
computational  methods  seems  quite  useful  and  promising  for  a  variety 
of  purposes .  First,  the  computational  methods  can  feasibly  be  used  to 
obtain  numerical  solutions  to  specific  problems  beyond  the  range  of 
existing  analytical  results.  Second,  the  exploration  of  numerical 
examples  may  provide  insights  concerning  problem  structure  that  could 
suggest  further  analytical  extensions.  Third,  computers  can  be  used  to 
do  sensitivity  tests  so  that  the  relative  advantages  of  alternative 
solutions  can  be  assessed;  thus,  an  approximate  solution  may  not  differ 
too  much  from  an  optimal  one,  but  it  may  be  much  simpler  to  obtain. 
Fourth,  by  Monte  Carlo  sampling,  the  computers  can  be  used  to  estimate 
variances  of  expected  value  solutions. 

This  hfemorandvun  deals  chiefly  with  the  first  of  the  purposes;  and, 
as  part  of  the  numerical  solutions,  it  considers  the  fourth  purpose 
as  well.  We  endeavored  to  compute  the  mean  values  of  selected  importsint 
inventory  characteristics:  number  of  shortages  and  number  of  overages 
per  time  period.  This  problem  can  be  solved  analytically  for  the  (S,s) 
Inventory  policy  model  with  a  zero  procurement  lead-time  and  exponential 
demand;  however,  for  the  non-zero  cases,  analytical  solutions  have  not 
been  derived. 

We  therefore  present  a  statistical  method  applicable  to  both  zero 
and  non-zero  procurement  lead-time  cases,  for  computing  the  expected  or 
mean  value  of  shortages  and  overages  per  period.  Then  we  apply  the 
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method  to  these  cases .  Our  intent  in  applying  this  statistical  procedure 
to  the  zero  case,  despite  our  ability  to  solve  it  analyiiically,  is  to 
check  whether  the  method  produces  estimates  of  the  true  mean  value  which 
have  some  precision  and  confidence.  We  check  this  approach  by  apply¬ 
ing  the  statistical  method  to  a  very  large  number  of  independent  trials 
and  make  an  estimate  of  mean  shortages  and  overages  each  time.  The  fre¬ 
quency  distribxxtion  of  these  estimates  can  then  be  used  to  evalmte  the 
extent  to  which  the  specified  precision  and  confidence  are  realized 
throiigh  use  of  the  statistical  method. 

We  have  also  used  a  similar  approach  with  the  non-zero  procurement 
lead-time  cases,  but  we  do  not  have  the  true  value  of  the  mean  ntunber 
of  shortages  and  overages.  As  a  substitute  for  these  true  values  we 
have  used  a  statistical  estimate  of  each  mean,  based  on  several  hundred 
observations .  Therefore,  the  evaluation  of  the  estimating  procedure 
is  not  as  rigorous  as  for  the  zero  procurement  lead-time  case;  the 
findings  obtained,  however,  should  give  an  empirical  measure  of  the 
validity  of  the  estimating  procedure  for  the  non-zero  case. 

The  study  thus  tests  and  evaluates  the  statistical  procedure  used 
to  estimate  mean  shortages  and  overages ,  It  is  a  companion  piece  to 
an  earlier  RAMD  Memorandum*,  subsequently  referred  to  as  RM-52if2,  which 
presented  the  method  for  computing  such  sanqile  sizes  . 


* 

M.  A.  Geisler,  The  Sizes  of  Simulation  Samples  Required  to  Compute 
Certain  Inventory  Characteristics  with  Stated  Precision  and  Confidence, 
The  RAND  Corporation,  RM-52J4-2-H?,  October,  1962. 
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II.  AMLYTIC  RESULTS  DERIVED  FOR  THE  ZERO  EROCUREMEMT  lEAD-TlME  CASE 

We  first  consider  an  inventory  model  with  a  zero  procurement  lead- 
time  governed  hy  (S,s)  policies.  We  assume  that  a  particular  set  of 
values  (S,s)  has  been  selected,  so  that  whenever  the  stock  level  x  falls 
to  or  below  s,  positive  ordering  is  immediately  enacted  to  raise  the 
level  to  S  upon  delivery.  t?hen  the  quantity  of  goods  in  supply  x 
exceeds  s,  then  no  ordering  is  done.  We  allow  x  to  assume  any  possible 
real  value.  A  negative  stock  level  should  be  interpreted  as  the  amount 
owed  to  consumption.  Thus,  all  demand  will  be  ultimately  satisfied, 
and  it  is  therefore  meaningful  to  refer  to  negative  stock  levels.  We 
also  assume  that  the  density  of  demand  f(§)  is  known,  so  that  in  each 
time  period,  a  demand  §  has  probability  f(5)  of  occurring.  Then,  if 
x^  =  stock  level  at  end  of  period  n,  we  have: 

x-§  -if  s<x£S 

n  *  n 

S  -  5  if  X  s  s  . 

n 


4 


Var  (x^)  = 


12^  3 


1  +  X  A 


T 


^  [x^A^  -  2XA  -  6 


X  X. 


iL  A^  +  i  a3  +  A®  +  ^ 

•\r\  ^  -x  ii  4i  ^ 
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A  + 


and 


P  =  P^ 

^n\+p  W+1 


vhere  A  =  S  -  s . 

These  results  are  derived  in  M.  A.  Geisler,  "Some  Statistical 
Properties  of  Selected  Inventory  Models, "  Naval  Research  Logistics 
Quarterly,  June,  1962.  The  article  also  derives  the  following  results 
of  interest  to  this  tfemorandxim.  If  y^  =  shortages  in  the  n-th  period. 


we  have: 


'n 


0  if  X  s  0 
n 


~x  if  X  <  0. 
n  n 


For  y^,  we  have  the  following  analytic  results: 


-Xs 


X(1  +  XA) 


9 


5 


)  . 


-Xs  -Xs 

Var  (jJ  =  - -  (2  -  j-TT-j^)  , 


“  X^(l  +  \h) 


Vn+l 


e~^^(l  ->-  XA)  ~ 

2(1  +  xa)  - 


Since  Is  not  a  Markov  process,  we  cannot  infer  p  ,  p  =  2,  5,  ••• 

^n^n+p 

from  p  ,  as  we  covild  with  x  , 

Vn+1  ^ 

The  term  "overage",  refers  to  the  positive  amount  of  stock  left 
at  the  end  of  the  period  before  ordering.  If  v^  =•  overage  In  n-th  period, 
then; 


I  x„,  -if  >  0 

I  0,  if  s  0, 


For  v^,  we  have  the  following  results:  true  for  each  n: 


”^+  X  /o2  2\ 

n^  X{l  +  Xa)  2(1  +  Xa) 


Var  (v  )  -  t.  -g  ■- +  UfL--5^. 

X^(l  +  XA)  5d  +  xa) 


~\s  +  -  1  ^  X 

.  x(i  +  u)  2(1  +  xa) 


(S^  -  s^)  I  , 


6 


-\s 


+  e 


-\s 


+  X  /S^  -  1  /S^  -  S^v 

I+U''  5  ^"1+\a'‘  2  ^ 


1 

1  +  \A 


rxs  +  -  1  ^  X 

L  X(1  +  XA)  2(1  +  XA) 


Finally, 


V  V  .  T 

n  n+1 


Var(vj^) 


The  expressions  given  above  for  Cov  Vax(v^)  can 

then  be  substituted  in  p  to  get  an  explicit  solution  for  p 

n^n+1  ^n\+l 

in  terms  of  X,  s,  and  A.  Here,  too,  v^  is  not  a  Markov  process  so  that 

ve  cannot  infer  the  behavior  of  p  ,  „  ,  from  p 

n  n+p  ^  n  n+1 

To  sxmnarize,  the  analytic  derivation  of  results  for  the  zero 

procurement  lead-time  case  with  exponential  demand  is  con^jlete  for 

x^  in  that  we  can  compute  exactly  the  mean,  variance,  covariances  for 

all  time  lags,  and  the  corresponding  correlation  coefficients.  For  y^ 

and  v^,  however,  we  can  compute  analytically  only  the  mean  and  variance, 

and  the  covariance  and  correlatiGn  coefficient  for  only  one  period  lag. 
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III.  AHALYTIC  RESUIffS  DERIVED  FOR  THE  NON-ZERO 
HtOCUREMENT  LEAD-TIME  CASE 

Analj^tically,  we  can  say  much  less  about  the  properties  of  the 
non-zero  procurement  lead-time  case.  We  have  been  able  to  do  the 
following.  If  we  let  z^  =  sum  of  on-order  plus  on-hand  stock  in  n-th 
period,  before  ordering,  we  then  can  show  that  z^  for  the  non-zero 
procurement  lead-time  case  is  equivalent  to  for  the  procurement  lead- 
time  case.  We  can  thus  compute  exactly,  for  z^,  its  mean,  variance, 
covariance  for  all  time,  lags,  and  the  corresponding  correlation  co¬ 
efficients  .  Also,  the  expressions  are  identical  with  those  given  above 

for  x^,  replacing  x^  by  z  . 
n  °  n  n 

We  can  ssj  further  that  if: 

x^  =  on-hand  stock  level  at  end  of  period  n  (which  can 
assinne  any  real  number  value), 
t  =  procurement  lead-time,  measured  in  number  of  time 
periods  from  order  to  delivery, 
then  the  following  relation  holds; 

*n+t-l  ~  ^n  ~  ^  ”  ' ' '  "  ^+t-2  ' 

where  5^^  =  demand  in  n-th  period. 

If  we  then  specialize  this  relation  to  the  case  of  t  =  2,  we 
obtain  Xq4.]_  “  2-period  procurement  lead-time  case, 

we  can  derive  the  following  results  analytically: 


f 


K(>„)  ■  E(z„)  -  ± 
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Vax  (x  )  =  Var  (z  )  +  , 

n  ^ 


2  2  2  2 
_  ,  X  „  .  /_  _  x  ,  \  ^  -  \S  -  \  s  +  3\S  -  5 

(“'n’ Vl^  ^  .  ^ 


“n’  n+1' 


XX,, 

n  n+1 


2  2  2  2 

„  /_  _  N  ^  X  3^=  -  XS  -  X  s  +  3XS  -  5 

Var  (z  )  + 

X 


These  results  are  meager  indeed,  so  ve  therefore  explored  the 
posslhillty  of  using  con^juting  techniques  to  augment  and  extend  them. 
In  the  following  sections,  we  present  the  statistical  technique  and 
the  resxxlts  obtained  with  it. 
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ry.  DESCRIETIOH  OF  STATISTICAL  -EECHNIQUE* 

We  would  like  to  compute  estimates  of  the  mean  nmber  of  shortages 
and  overages  occurring  in  an  inventory  model  operating  under  (S,s) 
policies .  We  consider  both  the  zero  and  the  non-zero  procurement  lag 
cases .  We  present  the  zero  lag  case  to  provide  an  absolute  assessment 
of  the  statistical  method  which  we  are  compelled  to  use  to  obtain  the 
properties  of  the  non-zero  lag  cases . 

ZERO  ffiOCUREMENT  lEAD-TBE  MOEEL 

We  are  interested  in  calculating  the  saiig)le  sizes  needed  to  assure 
that  the  sample  estimates  of  the  mean  shortages  and  overages  (per  time 
period)  differ  in  absolute  value  from  the  corresponding  true  value  no 
more  than  and  respectively,  with  95 "per -cent  confidence.  Thus, 
if  y  is  the  sample  estimate  of  the  mean  shortages  per  time  period  and  Y 
is  the  true  value,  then  we  want  to  find  that  minimum  value  of  n  ,  the 
san^ile  size,  such  that 

Er  {ly  -  Y|  >  kJ  £  O.O5  , 


Much  of  the  material  of  this  section  has  been  drawn  from  the 
companion  piece,  RM-52te,  to  provide  that  necessary  description  which 
will  help  the  reader  to  follow  the  test  procedure  and  results  presented 
later  in  this  Memorandum. 
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Similarly,  if  v  Is  the  sample  estimate  of  the  mean  overages  per 
time  period,  and  V  is  the  true  value,  then  we  want  to  find  that  minimum 
value  of  n^,  the  sarnie  size,  such  that 

Pr  ||v  -  V|  >  kJ  s  0.05  , 


-  n=l 

V  =  — - -  . 

V 

Subsequently,  where  ?  and  V  are  estimates  of  the  true  values  Y  and 
V,  we  let  =  ?  and  The  manner  of  estimating  ¥  and  ^  is 

described  in  Sec.  V.  Thus,  we  seek  those  minimum  saii5)le  sizes  for 
estimating  y  and  v  such  that  there  will  be  95-per-cent  confidence  that 
y  and  v  will  differ  from  their  corresponding  true  values  by  no  more 
than  approximately  the  tioie  value  in  absolute  amount.  This  is  the  same 
as  requiring  the  sample  sizes  to  differ  no  more  than  approximately  100 
per  cent  from  the  correspondingly  true  value,  with  95-per-cent  confidence. 

If  we  assume  y  and  v  have  approximately  normal  distributions, 
a  reasonable  assun^ition  in  view  of  the  Central  Limit  Theorem,  then  we 
can  use  Tschebycheff ’s  inequality,  and  choose  n  and  n  such  that 

V  • 

(1.96) cfy  s  Y  , 

(1.96) (t_  £  Y , 

for  attaining  a  ^-per-cent  confidence  level. 
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2  2 

We  can  now  calculate  cr_  and  ct_  using  the  following  formulas! 

y  V 


m 


aH  =  -^  {l  +  2  y  (l  -  P  )}, 
y  n  I  /L  '  m+l  ^p,y'P 


P=1 


m 


p=l 


where 

2 

a  =  variance  of  shortages  (per  period); 

•j 

2 

CT^  =  variance  of  overages  (per  period); 

P-,  „  =  P-th  order  lag  correlation  for  shortages; 

Jr  f  J 

p  =  p-th  order  lag  correlation  for  overages;  and 

ir9  ’ 

m  =  maximum  lag  for  which  correlations  are  computed;  so  that 
p  =  1^  2,  . . .  m. 


Hiese  expressions  are  modifications  of  the  more  usual  forms  for 
calcvilating  the  variance  of  the  mean  of  the  autocorrelated  series,  as 
given  by  Moran.*  They  allow  for  the  necessity  of  truncating  the  cal¬ 
culation  of  autocorrelations;  and  by  weighting  more  heavily  the  auto¬ 
correlations  for  the  shorter  lags,  they  attempt  to  reduce  the  likelihood 

2  2 

of  computing  negative  values  for  and  . 

2  2 

If  we  now  substitute  for  and  cr^  in  the  above  inequalities,  and 
solve  for  the  minimum  values  of  n  and  n  that  satisfy  these  Inequalities, 

y 


* 

P.  A.  P. 
York,  1959. 


Moran,  The  Theory  of  Storage, 


John  Wiley  and  Sons,  New 
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we  obtain: 


min(ny)  = 


Ij 


and 


rain(n^)  = 


where  [x]  =  greatest  integer  in  x. 

2  2 

As  stated  above,  we  can  compute  and  analytically  and  have 
therefore  used  their  analytical  values  in  computing  min  (n^)  and 
min  (n^).  We  cannot,  however,  compute  ^  and  ^  analytically. 
Instead,  we  have  to  estimate  them  with  Monte  Carlo  techniques,  thereby 
deriving  r  as  an  estimate  of  p  and  r  as  an  estimate  of  p 

^  p,y  p>y  p.v  ’^p,v 

The  technique  by  which  the  calculations  of  r  and  r  were  done  in 

!P^  y  p  ^  ^ 

a  particular  study  and  the  resulting  average  sample  sizes  with  their 
standard  deviations  are  presented  in  Tables  1  and  2  of  RM-52i+2. 


N0N-2ER0  roOCUREMENT  LEAD-TIME  MODEL 


The  approach  taken  to  conqjute  the  sample  sizes  for  this  model  is 
the  same  as  that  used  for  the  zero  jjrocurement  lead-time  model.  The 
main  problem  in  using  the  formulas  for  min  (n^.)  and  min  (n^)  is  that 


2  2 

we  cannot  derive  the  exact  values  0^  and  analytically.  Instead, 

2  2  2  2 

we  have  to  estimate  a  by  S  ,  and  a  by  S  .  The  procedure  for  doing 

y  y  V  V 

this  in  computing  sample  sizes  is  also  described  in  RM-32U2,  This 


procedure  was  applied  to  compute  sample  sizes  for  the  non-zero  procure¬ 
ment  cases  of  2-,  and  10-perlod  lead-times.  The  relevant  average 
sample  sizes  with  their  standard  deviations  are  given  in  RM-32i)-2:  Tables 
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5-6  for  the  2-period  procurement  lead-time,  Tables  9-10  for  the  5-period 
pi'ocurement  lead-time,  and  Tables  15-li+  for  the  10-period  procurement 


lead-time . 


V.  TEST  OF  SOIATISTICAL  OSCHNIQUE 


To  make  the  sample  size  calculations,  it  was  necessaiy,  in  effect, 
to  create  a  statistical,  model  of  the  sangiling  process .  This  was  pri¬ 
marily  represented  hy  the  formulas  specified  for  computing  min  (n  ) 

y* 

and  min  (n^) .  To  apply  these  formulas  we  had  to  make  three  decisions . 

First,  we  had  to  choose  an  appropriate  value  of  p,  the  truncation 

point;  and  second,  we  had  to  select  large  enough  Monte  Carlo  samples 

to  ensure  reasonably  valid  estimations  of  the  covariances  and  variances 

used  in  the  formulas.  Third,  although  theoretically  in  the  case  of 

exponential  demand  distributions  the  sample  size  computations  are 

Independent  of  the  value  of  x  used,  we  wished  to  use  one  value  of  \  in 

the  exponential  demand  distribution,  and  have  it  apply  to  the  range  of 

values  used  for  s  and  A.  We  chose  the  value  of  \  =  1.  Also,  we  were 

assuming  that  the  normal  distribution  assun^tion  could  be  used  in 

coir^juting  confidence  intervals  for  Y  and  V,  and  that  this  assumption 

would  also  hold  for  coii5)utlng  values  for  niin  (n  )  and  min  (n  )  that 

y  V 

would  result  in  the  specified  precisian  and  confidence  levels .  Thus, 
a  number  of  fairly  significant  but  subtle  assumptions  and  decisions 
were  required  to  establish  the  basis  for  utilizing  the  formulas  given 
above  for  min  (n^)  and  min  (n^). 

TEST  OF  ZERO  HtOCUREMEHT  I£AD-TIME  CASE 

To  evaluate  the  statistical  technique,  we  then  reversed  the  pro¬ 
cedure  used  in  RM-.5242  to  compute  min  (n  )  and  min  (n  ) ,  We  used  the 

J  * 

limiting  density  for  x  =  stock  level,  given  by: 
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cp(x)  = 


f  -  L_ 

1  +  XA 


Xe 


-X(s-x) 


1  +  XA 


if  s<x^A  +  s 


if  X  £  s. 


For  specific  values  of  X,  L,  and  s,  cp(x)  is  defined.  Using  this 
density  function,  a  value  of  x  was  selected  randomly,  and  set  equal 
to  Xq.  Then,  using  the  value  of  X  =  1,  a  specified  number  k  of  random 

y 

samplings  of  was  made  using  the  demand  density  function  f(g)  =  e”-, 
so  that  §Q,  ...  5-  were  obtained.  The  san5)le  size  k  was  selected 

Ky-t  y 

so  that  it  would  be  equal  to  the  mean  sample  size  as  given  in  Table  1 
of  RM-5242,  for  the  specified  s  and  A. 

Then,  using  the  transition  equations: 

'^n'^  if  s<x^^A  +  s 

^n+1  I 

^  if  Xn  s  s, 


we  calculated  the  random  variables  x^...  x^^  ,  thus  providing  a  random 

y 

sequence  of  k^.  values  for  x^  from  which  we  calculated 


0  if  X  >  0 
n 


yn  = 


-X  if  X  <  0. 
n  n 


Then  we  computed 


k 


it 


k  " 

y 


l6 


which  gave  the  mean  number  of  shortages  experienced  in  that  particular 
sample . 

We  repeated  the  process  to  compute  the  mean  number  of  overages,  v. 


From  Table  2,  we  obtained  the  requisite  mean  sample  size  k  .  If  k  2  k  , 

V  y  ^ 

vre  used  the  first  k^  random  samplings  of  ...  taken  from 


.  If  k  >  k  ,  however,  we  then  added  random  samplings 
u  K  "X  V  y 

y 

”k  ’  ■  ■  ‘  'k  1  earlier  series  ^  ,  using  the  same 

y  v“  y"-^ 

starting  value  Xq  drawn  from  cp(x) ,  We  could  thus  obtain  x^,  . . .  , 


from  which  v^  could  be  computed  by  means  of: 


if 


X  ^  0 
n 


v  = 
n 


if 


n  <  0. 


This  gave  v^,  . , 


Vj^  ,  from  which  we  confuted 
^  k 

V 


I 


_  n=l 
V  =  — r- 


V 


And  this,  in  turn,  gave  the  mean  number  of  overages,  v,  in  that  particular 
sample . 

To  test  the  precision  and  confidence  such  a  statistical  procedure 

would  give,  we  repeated  the  above  procedure  1000  times,  selecting  Xq  from 

cp(x)  and  a  sample  of  values  for  each  of  these  1000  random  samples . 

Vfe  thus  obtained  1000  values  of  j  and  v  for  the  specified  sample  sizes 

of  k  and  k  ,  respectively.  In  RM-5242,  k  and  k  were  computed  to 
y  V  y  V 

produce  sample  estimates  of  mean  shortages  and  overages  differing  no 
more  than  approximately  100  per  cent  from  the  corresponding  true  mean 
values,  with  95-per-cent  confidence.  ^Fbrther,  for  the  zero  procurement 
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lag  case,  we  know  the  true  values  of  the  mean  shortages  and  overages, 
which  we  designate  hy  Y  and  V  respectively,*  We  could  then  compute: 


for  each  of  the  1000  samples.  If  for  a  given  sample. 


then  we  had  an  estimate  which  satisfied  the  precision  condition  of 
being  within  100  per  cent  of  the  true  value .  We  could  then  count  the 
number  of  samples  out  of  the  1000  independent  ones,  say  N  ,  satisfying 

y 

this  condition.  Then, 

X  100  per  cent 

gave  the  confidence  realized  with  the  estimating  procedure.  This  same 
procedure  was  followed  with  the  mean  number  of  overages, 

N 

X  100  per  cent 

and  the  corresponding  confidence  values  computed. 

This  test  was  done  for  7  pairs  of  values  (s\,  L\)  that  seemed  to 
encompass  interesting  inventory  policy  cases .  The  results  of  this 
test  appear  in  Table  1  In  which  we  show  two  Interesting  statistics: 
first,  the  confidence  levels  actually  achieved  (with  95-per-cent  con¬ 
fidence  expected),  and,  second,  the  comparison  of  the  true  mean  values 
for  shortages  and  overages  with  the  estimate  obtained  from  averaging 
the  sang)ls  means  over  the  1000  sauries  computed.  To  be  more  precise, 

*These  values  of  Y  and  V  are  given  in  RM-5242,  Tables  5  and  ^ 
respectively. 
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if  and  are  the  mean  shortages  and  mean  overages  estimated  from 
the  i-th  sample,  where  i  =  1,  ...  1000,  then  we  confuted 


1000 

1000 

t  h 

and 

I 

i=l 

i=l 

1000 

1000 

Table 

1 

EVAUJATIOK  OP  STATISTICAL  ESTIMTING  mOCEDURES 
Itercentage  Confidence  Attained,  and  Comparison  of  True  and  Estimated 
hfean  Shortages  and  Overages:  Zero  Procurement  Lead-Time  Case 


Inventory 
Policies 
(sX,  AX) 

Per  cent 
Confidence 
Attained 

ffean 

Shortages 

Xfean 

Overages 

Sample 
Sizes  Used 

Short^ 

Ovg^ 

Est*^ 

True 

Est 

True 

Short 

Ovg 

0.01 

0.1 

97.6 

95.6 

1.05 

0.90 

0,006 

0,006 

5 

36 

0.1, 

0.1 

96.7 

95.7 

0.96 

0.82 

0,02 

0,02 

6 

19 

0.1, 

5.0 

95.0 

98.9 

0.15 

0.15 

2.9 

2.2 

32 

2 

1,0, 

0.01 

95.3 

99.1 

0.37 

0.56 

0,h5 

0.57 

16 

5 

1.0, 

1,0 

9^.5 

99.8 

0.18 

0.18 

1.58 

0.93 

30 

2 

1.0, 

5.0 

97.1 

99.9 

0.06 

0,06 

5.8 

3.0 

115 

1 

3.0, 

0.01 

97.2 

99.9 

0.05 

0,05 

i^,l 

2.1 

197 

1 

^Shortages 

Overages 

*^Estiinated 


and  these  values  are  the  estimated  mean  shortages  and  estimated  mean 
overages  given  in  Table  1. 

From  Table  1,  we  see  that  the  percentage  of  confidence  actually  at¬ 
tained  in  the  test  came  reasonably  close  to  the  desired  95”Per-cent 
value,  the  results  for  the  shortages  ranging  from  9^.5  to  97*6  per 
cent,  and  for  the  overages  ranging  from  95-6  to  99 >9  pen  cent.  Of  the 
7  values,  6  exceeded  95  per  cent  in  the  case  of  the  shortages j  all  7 
exceeded  that  figure  for  the  case  of  overages  .  A  partial  explanation 
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for  the  tendency  of  the  confidence  level  to  exceed  95  per  cent  is  that 

in  selecting  min  (n  )  and  min  (n  )  from  Tables  1  and  2  respectively 
y  V 

of  RM-32ii-2,  the  values  given  there  were  rounded  upward  in  all  instances. 
Thusj  if  n^  was  35*^  in  Table  2,  then  36  was  used.  The  effect  of  rounding 
would  be  greatest  with  the  smaller  sample  sizes,  and  may  help  explain 
the  relatively  higher  confidence  values  found  for  the  overage  cases 
when  samples  of  2  and  1  were  used. 

Relatively  good  agreement  was  also  found  between  the  true  measa 
shortage  values  Y  and  those  estimated  from  the  1000  sample  values  of  y. 

We  note  that  the  agreement  was  better,  the  larger  the  sample  size  used 
to  estimate  y.  A  similar  finding  occurred  in  comparing  the  true  value 
of  mean  overages  V  and  the  mean  overage  value  estimated  from  the  1000 
sanqple  values  of  v. 

Thus,  from  the  data  contained  in  Table  1,  we  would  infer  that  the 
statistical  procedure  tends  to  do  what  it  was  devised  to  accomplish: 
it  produces  estimates  of  the  mean  value  of  shortages  and  overages  within 
100  per  cent  of  the  true  value,  with  95-pex-cent  confidence. 

TEST  OF  NON-ZERO  IROCUREMEMT  lEAD-TIME  CASES 

A  procedure  similar  to  that  described  for  the  zero  procurement 
lead-time  case  was  followed  for  the  non-zero  procurement  lag  models  of 
2,  5^  and  10  time  periods.  As  mentioned  in  Sec.  Ill,  if  z^  stands  for 
the  on-hand  plus  on-order  stock  in  the  n-th  period  (before  ordering), 
then  z^  has  the  same  statistical  properties  in  the  non-zero  procurement 
lag  case  as  x^  in  the  zero  lag  case.  Thus,  for  the  exponential  distri¬ 
bution,  f(5)  =  \e"^^: 
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cp(z)  = 


1+XA 

^^-X(s-z) 

1+\A 


if  S<Z:SS  +  A 


if  Z  £  s . 


In  addition,  if  x  =  the  amount  of  stock  on  hand  at  start  of 
n 

period  n,  including  stock  delivered  in  period  n,  and 
=  demand  in  n-th  period,  and 
t  =  procurement  lead-time  in  number  of  time  periods  from 
order  to  delivery,  then,  as  stated  in  Sec.  Ill: 


:  , ,  ,  =  z  -  ?  “  •  •  •  ~  S  .  ^  for  t  5:  2 . 

n+t-1  n  ^n  ^n+t-2 


We  thus  could  select  Zq  from  c)(z),  then  select  the  appropriate 

series  of  §  values,  and  derive  the  series:  y^ ,  . . .  y,  for  shortages 
n  j.  K 

y 

and  v^,  . . ,  for  overages,  paralleling  the  procedure  described  in 

V 

this  section  for  the  zero  lag  case.  We  thus  computed 


y  = 


y 


and 


V  = 


k 

V 


I 


V 

n 


n“l 


for  the  particular  saii5)le  of  values  Zq,  5^,  ...  ^  where  k  =  k^  or 

k  =  k^  depending  on  which  is  larger.  We  then  repeated  this  sampling 
and  computing  process  1000  times,  obtaining  1000  independent  values 
each  of  y  and  of  V.  In  the  non- zero  lag  cases,  we  did  not  know  the  true 
values  of  Y  and  V,  the  mean  values  of  shortages  and  mean  overages  per 
period,  respectively.  Instead,  we  used  approximate  values  designated 
by  If  and  V,  which  were  computed  from  a  time  series  of  5OO  vaulues  of 
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and  v^.  The  values  of  Y  and  used  are  given  in  RM-52^2:  Tables 
7  and  8  for  the  2-period  lag  case,  Tables  11  and  12  for  the  5"Period 
lag  case,  and  Tables  15  and  l6  for  the  10-period  lag  case.  With  this 
information  we  could  compute 


Y  Y 


for  each  of  the  1000  sauries,  and  establish  whether 

!^|  ^  1  and  l^i  ^  1 
Y  V 

in  each  sample .  We  could  then  compute 

N  N 

yr  -y 

lote  1000 

which  gave  the  percentage  of  samples  for 

tv  -  Yt 

s  1  and  -  ^  1  ^ 

Y 

respectively.  These  two  percentages  thus  gave  an  indication  of  the 
confidence,  within  100  per  cent  of  the  approximate  true  values,  that 
the  statistical  method  would  yield  in  its  estimates  of  mean  shortages 
and  overages ,  We  also  coioputed  the  average  of  the  y  and  of  the  v  values 
for  the  1000  samples  to  determine  how  well  these  averages  agreed  with 
the  corresponding  approximate  true  value  Y  and  Y. 

The  results  of  these  calculations  for  selected  inventory  policies 
are  given  in  Table  2  for  the  2-period  lead-time  case,  in  Table  3  for 
the  5-period  lead-time  case,  and  in  Table  for  the  10-period  lead-time 


1^1 


case , 
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The  results  found  in  each  of  the  non-zero  lead-time  cases  seem 
to  a^ee  well  with  those  for  the  zero  lead-time  case.  That  is^  the 
confidence  actxially  attained  with  the  comfiuted  sample  size  is  not  far 
from  the  expected  value  of  95  pei*  cent^  with  most  of  the  confidence 
values  exceeding  that  figure.  The  highest  confidence  values  tend  to 


Table  2 

EVALUATION  OF  STATISTICAL  ESTIMATING  IROCEDURES 
Percentage  Confidence  Attained,  and  Cougar ison  of  "Approximately  True"  and 
Estimated  Lfean  Shortages  and  Overages:  2-rteriod  Procurement  Lead-Time  Case 


Inventory 
Policies 
(s\,  dX) 

Per  Cent 
Confidence 
Attained 

Mean 

Shortages 

Mean 

Overages 

Sample 
Sizes  Used 

Short^ 

Ovg^ 

Est*^ 

Approx 

True<i 

Est 

Approx 

True 

Short 

Ovg 

0.01,  0.01 

97  A 

e 

2.0 

2.0 

n 

5 

0.01,  J+.o 

97.3 

CO 

0.7 

0.6 

1.0 

20 

7 

1.0,  0.01 

97.2 

99.9 

1.2 

1,1 

0,1 

0.2 

13 

32 

1.0,  0.1 

98.5 

1.2 

1.0 

.. 

.. 

I4 

2.0,  0.01 

98.9 

99.9 

0.7 

0.5 

0.7 

0.5 

29 

12 

2.0,  3.0 

— 

99.9 

— 

2.0 

2.1 

3 

3.0,  10.0 

-- 

99.9 

— 

— 

7.0 

6.4 

— 

2 

s  (i 

Shortages  Approximately  true 

b  e 

Overages  Not  calculated 

‘Estimated 


be  found  with  the  smallest  sample  sizes;  this  occurrence  is  believed 
to  be  attributable  to  the  author’s  practice  of  reminding  the  average 
sample  sizes  in  RM-5242  upward  to  the  next  integer  value.  Also,  the 
estimated  mean  shortages  and  mean  overages  conpited  from  the  1000  samples 
of  specified  size  seem  to  agree  quite  well  with  the  approximately  true 
mean  shortages  and  overages,  especially  for  the  larger  specified  sainple 


sizes . 
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Table  3 


EVALUATION  OF  STATISTICAL  ESTIMTING  TROCEDURES 
Percentage  Confidence  attained,  and  Comparison  of  "Approximately  True "  and 
Estimated  tfean  Shortages  and  Overages:  5“P2riod  Procurement  Lead-Time  Case 


Inventory 
Policies 
(sX,  A\) 

Par  cent 
Confidence 
Attained 

Mean 

Shortage  s 

Mean 

Overages 

Samp 

Sizes 

le 

Used 

Short®' 

Ovg^ 

Est° 

Approx 

True^ 

Est 

Approx 

True 

Short 

Ovg 

0.01,  0.01 

98.3 

e 

4.9 

5.0 

5 

... 

0.01,  3.0 

96.2 

— 

2.6 

3.2 

— 

-- 

11 

— 

0.01,  50.0 

.. 

99.9 

20.0 

20.9 

— 

6 

1.0,  2.0 

96.6 

— 

2.8 

2.8 

— 

— 

13 

— 

2.0,  0.01 

98.0 

— 

3.0 

3.0 

13 

2.0,  2.0 

97.7 

— 

1.9 

1.9 

— 

— 

22 

— 

4.0,  3.0 

-- 

99-4 

— 

-- 

1.4 

1.5 

— 

20 

4.0,  5.0 

— 

99.9 

— 

— 

3.0 

2.3 

11 

5.0,  4.0 

— 

99.8 

— 

-- 

2.9 

2.7 

10 

10.0,  0.01 

— 

99.9 

— 

-- 

4.0 

5.0 

IBS 

5 

£L  d 

.Shortages  ^ Approximately  true 

Overages  ®Not  Calculated 

'^Estimated 


Table  4 

EVALUATIOCr  CF  STATISTICAL  ESTIMATING  mOCEDURES 
Psrcentage  Confidence  attained,  and  Comparison  of  "Approximately  True"  and 
Estlnated  Mean  Shortages  and  Overages: 10-Period  Procurement  I^ad-Time  Case 


Inventory 

Policies 

(sX,A\) 

'  '  1 

Bar  cent 
Confidence 
Attained 

Mean 

Shortages 

Mean 

Overages 

Sample 

Sizes  Used 

Ovg^ 

Est*^ 

Approx 

True^ 

Est 

Approx 

True 

Short 

Ovg 

0,01,  0.01 

99.7 

10.0 

10.0 

5 

„ 

0.01,  5.0 

96.0 

5.2 

7.1 

16.5 

16.7 

9 

9 

1.0,  3,0 

98.6 

7.2 

6.9 

— 

10 

— 

3.0,  4.0 

98.5 

4.7 

4.6 

— 

— 

18 

— 

4.0,  0.01 

98.3 

5.9 

6.0 

— 

— 

13 

— 

4,0,  50.0 

— 

95 

— 

— 

19.6 

20.0 

— 

5 

5.0,  50.0 

99.9 

— 

— 

22.2 

25-6 

— 

6 

10.0,  10.0 

99.9 

— 

-- 

8.3 

10.4 

— 

U 

50,0,  10.0 

__ 

99.9 

— 

— 

45.3 

45.4 

-- 

1 

d 

^Shortages  Approximately  true 

^Overages  Calculated 

'Estimated 


2h 


Thus,  on  the  basis  of  this  test,  we  can  further  state  that  the 
estimating  procedvire  presented  in  RM-3242  for  conqcjuting  sample  sizes 
of  the  (S,s)-type  inventory  models  seems  to  be  valid,  in  that  it  does 
provide  the  confidence  desired  in  the  estimates  of  mean  shortages  and 
mean  overages. 


